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ABSTRACT 

A systematic technique is presented for modeling crack growth load 
interaction effects due to spectrum loading. The Multi-Parameter Yield Zone 
(MPYZ) model accounts for crack growth retardation, acceleration, and under- 
load effects. The load interactions are attributed to the residual stress 
intensity due to the plastic deformation at the crack tip. As part of an 
ASTM E24. 06.01 round-robin effort, fatigue crack growth was predicted and 
compared with test data for a variety of spectrum loadings. 
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effective applied minimum stress intensity after load 
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overload maximum stress intensity that is still effective prior to 
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the minimum stress intensity associated with the K cycle, 
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the residual stress intensity associated with the MPYZ load 
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interaction model, MN-m 

the residual stress intensity associated with the Willenborg 
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retardation model, MN-m 
current underload stress intensity, MN-m 
material crack growth constant 


stress ratio, S . /S 
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eff eff 

stress ratio after load interaction, K . /K 

mm max 

ratio of current underload stress to maximum overload stress, 
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UL max 

stress state dependent maximum * value of R 
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Aa 
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4k th 
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TT 
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maximum overload stress, MPa 

minimum stress associated with overload cycle, MPa 
underload stress, MPa 

value of 0 bejow which overload is unaffected 
value of 0 above which there is no overload effect 
plastic zone diameter associated with currently applied load, m 
plastic zone created by overload, m 


1 for plane stress 

3 for plane strain 


underload factor, ( Kpr - - K^) 

crack growth increment since overload applied, a QL - a ap' m 


stress intensity range, K - K . , MN-m 

^ max mm 


-3/2 


-3/2 


threshold stress intensity range at R = 0.0, MN-m 
threshold stress intensity range for a given R e ^, MN-m 
constant 3.141592654 
remote field stress, MPa 
material yield stress, MPa 


INTRODUCTION 

The early 1970's saw three significant concepts generated to deal with 
spectrum effects on crack propagation. They analytically accounted for the 
experimentally observed phenomenon of crack growth retardation [1-3] . These 
three approaches, commonly referred to as the Wheeler model, the Willenborg 
model, and the Elber (or crack closure) model, along with the creation of 
large and efficient digital computers, spurred the aircraft industry to 
analytically predict crack growth from realistic spectrum loadings. The 
Wheeler model is characterized by a shape-fitting parameter, "m" , that must 
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be experimentally determined for each condition of material, loading, and 
environment that is of interest. The Willenborg model, on the cither hand, 
does not incorporate an empirical shape factor but uses the material yield 
stress to give a plastic zone size, which is then used to calculate crack 
retardation. Although the Willenborg model does not incorporate a "fudge" 
factor, its predictions are not always satisfactory. The Elber model is 
based on the crack-closure phenomenon and uses an effective stress intensity 
which is based upon the stress needed to open the crack. The closure 
factor in the Elber model is determined experimentally from constant- 

^plitude data. The Elber model has been used in the aerospace community 
with limited success. These models have been modified by various investi- 
gators during the past five years with some degree of success [4-7] . 

Recent testing at General Dynamics [8-9] , as well as work performed by 
Hillberry [10-12] and others .[13-14] , have shown that various loading param- 
eters that are not included in the Wheeler or Willenborg models (such as crack 
growth acceleration and underloads) have a significant effect on spectrum 
crack growth. 

This paper will present a phenomenological load interaction model 
referred to as the Multi-Parameter Yield Zone (MPYZ) load interaction model. 

It is a subprogram within the computerized crack growth program CGR in ref- 
erences 15 and 16. The crack growth model is discussed first, followed by 
the load interaction model. The load interaction model is described in 
three sections: retardation, acceleration, and underload effects. The crack 

growth predictions are discussed and compared with the experimental results. 
The present' study was part of a round- robin effort that was sponsored by 
ASTM Subcommittee E24.06. The purpose of this round-robin effort was to 
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compare various methods for predicting crack growth through center cracks in 
2219-T851 aluminum alloy center-crack tension specimens subjected to aircraft- 
type spectrum loadings. The data used for the round robin were generated by 
Chang and Stolpestad [17] under an Air Force sponsored contract. All 
material data in this paper can be found in reference 17 , except as noted. 

CRACK GROWTH MODEL 

Numerous fatigue crack propagation models have been proposed in the 
literature to describe the relationship between the crack growth rate, da/dN, 
and the stress intensity range, AK. Many of these relationships consider 
such parameter as stress ratio, R, and fracture toughness, K c . Gallagher 
[Gallagher, J. P., "Fatigue Crack Growth Rate Laws Accounting for Stress 
Ratio Effects," ASTM Task Force E24.04.04, Report No. 1, 1974.] has summarized 
a number of these models. 

The Forman [24] crack growth model provides a good fit for the constant 
amplitude 2219-T851 aluminum crack growth data furnished for the spectrum 
crack growth prediction round- robin. For the present study, the Forman 
equation was slightly altered to the following form: 

da/dN = , 1 , 

(1 - R e£f ) K c - dK 

where 

m = 1 at R - 0 
m = 2 at R < 0 

The equation is exactly that suggested by Forman except for the 
m exponent applied to the stress ratio. The exponent m = 2.0 for nega- 

ef f 

tive R values. Equation (1) accounts for the crack growth behavior 
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by using the full stress intensity range, AK, and the corresponding R 

values, i.e., the stress range includes compressive portions. Load inter- 

ef f 

action effects are accounted for by adjusting R during each cycle, so 

e f f 

R does not necessarily equal the actual stress ratio of the current 

cycle. Most of the present paper is about the procedure used to find the 

ef f 

proper values of R to account for retardation, acceleration, and under- 
load effects. The parameters C and n were found from the supplied 
constant amplitude crack growth data [17] by a least-squares regression 
analysis. The critical stress intensity factor, K^, was chosen to be 
88 MN-m (80 ksi \Jin . ) , which resulted in C = 3.68 x 10 and n = 3.232. 

The author has formulated an opinion from observed crack growth 
behavior, and from crack closure considerations of Newman [19] , that the 
maximum effective stress ratio is a function of the constraint at the crack 
tip due to the different stress states, i.e., plane-stress or plane-strain 
or in between. Therefore a simple relationship has been used to determine 
the maximum allowable stress ratio used in equation (1) : 


R max = [ Z/t * 0-2]+ 0.6 (2) 


air \ a. 


is the plastic zone diameter for the applied K , 

max 


t is the material thickness, and a is 1.0 or 3.0 for plane stress or 


strain assumption, respectively. Notice that for plane-strain (Z/t + 0) 


so that R max = 0.6. Therefore if an applied cycle had a stress ratio equal 
to 0.75 9 hut vas cycled in a plane-strain region, a value of approximately 
0.6 would he used in equation (l) for R • . If the plastic zone diameter 
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equals or exceeds the thickness, then a state of plane-stress is assumed, and 

eff 

R = 0.8, i.e., maximum value for R in equation (1) is 0.8. 
max 

Figure 1 depicts the fit of equation (1) to the data. The solid line 
represents perfect correlation while the dashed lines represent a factor of 
two discrepancy. The constant amplitude data consisted of R values between 
-1.0 and 0.7 and maximum stress levels from 55.2 to 276 MPa. Figure 1 is a 
plot of the predicted crack growth rates versus the actual growth data used 
to evaluate the C and n in equation (1). Several constant amplitude 
tests at 276 MPa (40 ksi) were not included because of net section yielding; 
the author did not want to bias the crack growth equation for this extreme 
case. It is interesting to note that equation (1) does a good job of 
correlating the negative and positive stress ratios. 

Typically, many cycles in a variable-amplitude loading program produce 
stress intensity ranges, AK!s, that are so small no crack growth results. 

The threshold stress intensity range, Ak£^, is stress ratio dependent and 
is not accounted for in the modified Forman equation. If the applied Ak 
is less than Ak£ h , then equation (1) is not applicable and no crack growth 
results. The threshold for each cycle, Ak£^, was calculated from a read in 
value of Ak^ h at R = 0.0, as defined from constant amplitude data, such 
that 

AK TH = (1 " R6ff)AK TH (3) 

-3/2 

where AK^ is, from basic data [20], taken to be 3.3 MN-m for 
2219-T851 aluminum alloy. 

A method for assessing constant amplitude crack growth data as well as 
an efficient and accurate crack growth accumulation/integration routine for 
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the Forman equation have been previously presented by the author in ref- 
erence 21 and will not be repeated herein. 

In summary, then, to calculate an increment of crack growth equation (1) 
is used. Two situations may exist; AK is below threshold, and da/dN = 0; 

or AK is above threshold, and load interaction effects are taken into 

eff 

account by properly adjusting the effective stress ratio, R . Most of 

the remainder of this paper describes the procedures used to determine the 

eff 

proper values of R to be used in equation (1) during each load cycle. 

LOAD INTERACTION MODEL 

Realistic aircraft loading results in a mixture of high and low stresses 
associated with varying stress ranges. The influence of any given load cycle 
on crack growth behavior depends upon the prior load history, thus the term 
"load interaction" is used to describe these effects. The load interaction 
model used by the author to predict crack growth under spectrum loading is 
referred to as the Multi-Parameter Yield Zone (MPYZ) model. The MPYZ is a 
phenomenological model that is intended to account for such recognized load 
interactions as retardation caused by previous overloads, acceleration due 
to current overloads, and underload effects. Retardation results in less 
crack growth for a given load cycle than would be produced under the same 
constant amplitude condition. Conversely, acceleration results in more crack 
growth than would be expected under constant amplitude. The underloads 
result in a reduction of the retardation effects of the current overload. 

The MPYZ model began in 1973 [15] as a form of the Willenborg retarda- 
tion model [2]. During 1976, in reference 21, two parameters were added to 
the Willenborg model. These parameters correspond to ratios of a single 


8 



overload stress intensity to the following applied stress intensity that 
resulted in (1) no crack growth retardation (parameter A) or (2) complete 
crack growth retardation (parameter B) . Parameter B is very similar to the 
"shut-off" ratio implemented by Gallagher (16] . 

In the present model, load interaction effects are calculated utilizing 
a residual stress intensity, Kj^, concept- Crack growth retardation and 
acceleration are accounted for by decreasing or increasing, respectively, 

the effective stress ratio used in a modified Forman crack growth equation, 
ef f 

The R is a function of the residual stress intensity, k r , and the 

currently applied stress intensity. The effective stress range that is used 

to calculate crack growth remains the same, but the mean of the stress range 

is altered by to account for the appropriate load interaction. The 

eff 

stress ratio corresponding to a load interaction will be designated R 
and defined as follows: 


R eff 


K , 
min 


K 

max 


- Kr 


K eff 

min 

, eff 

max 


(4) 


where is the residual stress intensity which will be discussed in the 

following section. The will be positive or negative depending on whether 

the load interaction produces retardation or acceleration, respectively. The 

eff 

maximum allowable value of R is R as defined in equation (2) . A 

max 

simplified flow chart of the MPYZ model in the CGR computer program is 
presented in figure 2. The three primary phenomenological events of load 
interaction, i.e., retardation, acceleration, and underload effects, will be 
addressed in the following sections. 
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Retardation 


The amount of crack growth retardation for a given applied cycle depends 
upon the applied loading as well as the extent of crack growth into the over- 
load plastic zone, as suggested by Willenborg, Engle, and Wood [2] . 

Gallagher [5] expressed the Willenborg model as 



OL 

K 

max 



K max 


(5) 


W ^ OL 

where K is the Willenborg residual stress intensity factor; K is 

R ^ max 

the maximum stress intensity for the overload cycle; K max the maximum 

stress intensity for the current load cycle; Aa is the crack growth incre- 
ment between the overload cycle and the current cycle as depicted in fig- 
ure 3. The Willenborg model applies Z as the plastic zone radius 

associated with . But in the present MPYZ model, is defined as 

max OL 

the zone diameter such that 


Z 0L OtTT 


k ol 

max 


YS 


(6) 


where 

a - plastic zone constraint factor 
=3, plane strain 
= 1, plane strain 
O = material yield stress 

X b 

w 

The plastic zone diameter is used to calculate K^. Mills, Hertzberg, and 
Roberts [18] have shown that the retardation effects last approximately one 
plastic zone diameter in A514F steel. Himmelein and Hillberry [11] present 
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data that also imply the overload is effective over the plastic zone diameter 
for 2024-T6 aluminum. 

W 

The residual stress intensity, K , calculated in equation (5) physically 

R 

represents the difference between the stress intensity required to produce a 

current plastic zone equal to Z QL - Aa and the current applied stress 

w 

intensity , K . Retardation occurs when K is greater than zero. In 
J max R 

such cases the residual stress intensity at the crack tip is greater than 
that for an equivalent constant amplitude case, so there is less crack growth, 
i .e . , retardation. 

Figure 4 presents a typical aircraft load sequence. Notice that the 

overload and underload are labeled. The dashed line represents the decay in 

eff eff 

effective overload stress intensity, K . The value of K__ decreases 

OL Ob 

with each cycle of crack growth. 

Gallagher and Hughes [5] introduced a proportionality factor, <f) , 

R 

such that 

*R - vS < 7 > 

The term adjusts the amount of retardation according to several material 

and -load history parameters. Notice that the higher the value of <b the 

R 

more the retardation. 

The MPYZ load interaction model incorporates 

* R = C 1 - 0 - ‘ww] / [w - i-o) x d * 0 - v] (8) 


ii 


where 


^ = B/A; A = the ratio of K OL to K below which no retardation 

max max 

takes place (no-retardation ratio) 

B = the ratio of yP^ to K above which no crack growth 

max max 

is produced (shut-off overload ratio) 

R t = S riT /s2t,; ratio of current underload stress to current overload 
UxjI .nax 

stress. 

The [l.O - (Krn'^max^ factor accounts for the threshold level of 

crack growth. The (\p - 1.0) factor adjusts the retardation between the 

shut-off overload ratio B and the no-retardation ratio A as shown in figure 5. 

The (1.0 - R^) expression adjusts the amount of retardation according to 

the cyclic range of the spectrum. Similar spectra that differ only in the 

minimum underload will have different overall retardation. This phenomenon 

is also observed in era ;k ' iosure [19] where the maximum and minimum stresses 

in a spectrum usually set the crack opening stresses therefore controlling 

retardation. Notice that a spectrum with negative minimum (i.e., negative 

value of R t ) will produce less retardation for a given ratio of K° L /K 

^ " max/ max 

than a similar loading sequence with a positive minimum stress. 

Parameters A and B are considered to be material dependent param- 
eters. Ideally, A and B can be determined for a given material by con- 
ducting a series of well controlled systematic single overload tests. The 
ratio of the applied single overload stress to the following stress is varied 
from test to test in order to ascertain the limiting ratios that (1) produce 
no overload efiozz (i.e., ratio = A) and (2) produce no crack growth other 
than that produced by the overload (i.e.; ratio = B) . The types of tests con- 
ducted by Probst and Hillberry [12] are suited to this purpose. If these 
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systematic test results are not available for evaluation for a particular 
material, one has to try various combinations of A and B to predict 
crack growth and compare the predictions with available test results. Typical 
values of A and B found in the literature range from 1.0 to 1.5 and 
1.8 to 2.5, respectively, for various materials. 

The MPYZ load interaction model uses a simple algorithm to account for 
multiple overload effects [13]. Simply stated, the more overload cycles one 
has in succession the greater the retardation effect, i.e., the lower the 
value of B' . Thus this relationship is approximated for 2219-T851 aluminum 
by 

B ' ■ hvi + 2 -° 

where 

B' = the value of B to be used in equation (8) 

N ql = the number of overloads applied in succession 

The 2.0 is considered to be a lower bound for B' after many overloads 
for 2219 aluminum alloy 

Although the equation is simple, it accounted for the multiple overload 
effects in reference 17 quite well. 

Acceleration 

The crack growth during an overload cycle has been observed to be larger 
than one might expect from constant amplitude data [14] . This phenomenon is 
referred to as crack growth acceleration. The acceleration may be attributed 
to the crack and the resulting crack tip plastic zone growing into an area 
with less residual stress intensity than that of an equivalent constant 
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amplitude case, so that crack growth is less inhibited. For example, imagine 

that the plastic zone diameter, Z , in figure 2 was caused by a large 

ap 

enough load such that 


Aa + Z 

ap 


> Z 


OL 


( 10 ) 


This would result in a crack growth acceleration for the applied cycle. 

Equation (4) is also utilized for acceleration load interaction. When 

W 

inequality (10) is true, acceleration occurs, and K r is a negative value. 

The MPYZ load interaction model incorporates an acceleration adjustment 
W 

term, 4> A , for K r defined as 

4> a = (1.0 - <U) 

where the term (1.0 - R^) adjusts the amount of acceleration depending on 
the ratio of the current underload stress to the overload stress. Thus the 
residual stress intensity used in equation (4) is defined as 


K 


R 


y A R 


( 12 ) 


for acceleration. Note that <J>^ is used to calculate when is a 

eff 

negative value. Thus, the effective stress ratio, R , will increase. 


Underload Effects 

Occasionally in spectrum aircraft-type loading, a rather significant 
compressive load or a tensile load occurs that is lower than most of the 
previous minimum loads. This load, S^, is defined as an "underload," as 
depicted in figure 4. An underload has dramatic influence because it can 
significantly reduce retardation. The author has conducted tests to establish 
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the effect on crack growth caused by deleting the compressive stresses from a 
fighter cycle-by-cycle spectrum. These tests were conducted on 7475-T7651 
aluminum in Sump Tank Water (STW) at a maximum spectrum stress of 214 MPa 
(31 ksi) . The fighter spectrum contains relatively few compressive loads for 
the wing root bending moment spectrum. Those that are present have magni- 
tudes less than 15 percent of the maximum tensile load. As shown in figure 6, 
the deletion of these few small compressive stresses resulted in a 50 percent 
increase in life. This substantial increase in crack growth rate cannot be 
attributed solely to the compressive load excursion increasing the stress 
range and producing a negative stress ratio. Rather , the negative stresses 
tend to negate the crack growth retardation caused by previous tensile 
overloads . 

Hillberry and associates [10-12] have conducted several studies into this 
particular phenomenon. In reference 10 underloads of varying amplitude were 
applied iirenediately after an overload to assess the effect of each underload 
amplitude on the following crack propagation behavior. These data are very 
useful in modeling the underload phenomenon. Chang and Stolpestad [17] also 
have conducted several simple variable amplitude tests that lend themselves 
to characterizing the effects of underloads. 

In actual aircraft load histories, a significant underload may not 
immediately follow a large overload. Therefore, in order to generalize the 
model to realistic spectrum loadings, an effective overload is computed when- 
ever an underload cycle is applied. The effective overload represents the 

remaining contribution of K before underload application, to simulate 

max 

the conditions of Hillsberry's data upon- which the following model is based. 
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The effective overload is assumed to be the load required to create a plastic 


zone size of Z - Aa in figure 3. The effective overload stress intensity 
OL 

is therefore equal to where 

=K 0L /l-§M 1/2 (13) 

eff max ^ Z QL J 

The underload is accounted for by reducing the effective overload stress 
intensity factor using the ratio 


= V ~ V (141 

OL 

K eff - ^JL 

where K , K eff' and ^UL are shown in fi 9 ure K p r is tiie minimum value 

of stress intensity before the overload. Thus a new effective overload stress 

intensity factor , denoted K ' , is computed and used in accounting for crack 

OL 

growth retardation. The value of this effective overload stress intensity, 

K* , is expressed as 



K °ff 

Z“Y (Z " B) + K max 


(15) 


where Z is the value of 3 above which k’ is K and Y is the 

OL max 

value of 3 below which k'_ = K°^, . Thus, Y determines when the underload 

ol err 

does not reduce the amount of retardation, and Z determines when no retarda- 
tion occurs. The relationship among these variables is indicated in 
figure 7. 

Notice, in figure 7, that when 8 <• Y the underload has no effect on 
retardation and equals . Likewise, when 6 is greater 
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than Z, K„ equals K so no retardation results. If 8 is between 

OL nax 

Gf f 

the values of Y and Z the effective overload, K , is reduced. The 

OL 

relationship between B and depends on the values of Y and Z and 

can be represented by a straight line. 

The Y and 2 for the particular material must be determined in a 
similar fashion as that for A and B. Carefully controlled crack growth 
tests can be run in a systematic manner to completely determine the crack 
growth load interaction behavior. Tests as conducted in reference 10 are 
ideal for the determination of Y and Z. 


RESULTS OF SPECTRUM CRACK GROWTH PREDICTIONS 
The crack growth prediction computer program, CGR-LaR C (an updated 
version of CGR [15,15]), was used to make predictions for the present round- 
robin tests. The program required the usual descriptions of shape and 
material properties of the test specimen to be analyzed. In addition the 
program required the four load interaction parameters A, B, Y, and Z as 
discussed in the previous section. These parameters, as well as the material 
data for the 2219-T851 aluminum, are presented in table 1. Notice that 
A = 1.0 and Y = 0.0 such that there are essentially only two load inter- 
action parameters used to describe 2219-T851 aluminum, B and Z. The 
B = 2.3 agrees with data presented in reference 12. 

Reference 17 contained many simple single overload tests. The tests 
were not systematic (i.e., each test conducted at an incrementally different 
value of overload ratio) but they did cover a wide range of overload ratios. 


K 


.OL 


max 


K 


The A, B, Y, and Z terms were selected from the simple overload 


max 
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test observations and to provide the best fit of the predicted to the experi- 

mental test results for these simple tests. 

The round-robin program consisted of 13 tests. Five basic aircraft type 
spectra were used. Thr6e of the spectra were applied at three different load 
factors (i.e., different scaling of the stresses) and two spectra were applied 
at two different load factors. There was only one test per individual 
spectrum. More details are available in reference 17. 

The round-robin specimens contained a through center crack. The secant 
correction factor was used to calculate the stress intensity factor [23] : 



where 

0 = the remote stress 
a = half surface crack length 
b = half specimen width. 

Figure 8 is a plot of the predicted lives versus the test lives from the 
13 tests in the round-robin effort. The solid line represents exact predic- 
tions. The dashed lines are a factor of two high or low. The mean and 
standard deviations of the ratio of predicted life to test life are 0.97 and 
0.24, respectively. The standard error for the 13 test predictions is 


2 


N , 
_ pred 

^test 


S.E. 


13 


0.18 



Figure 9 presents the predicted and experimental crack growth curves of 
specimens M-90, M-91, and M-92 for comparison. The shape of the predicted 
crack growth curve is in excellent agreement with that of the experimental. 

All of the predictions were within a factor of 1.5 of the test results. 

In general, all predictions were extremely good with M-88 and M-89 being the 
poorest. (See table 2.) Notice in this table that the ratio of the predicted 
life to the test life is shorter at higher load factors within each spectrum 
(81-83, 84-86, 88-89, 90-92, 93-94). This is attributed to the changes in the 
constraint at the crack tip. All of the presented predictions assumed plane 
strain conditions in order to be slightly conservative overall. The higher 
the stresses the more likely plane stress occurs, especially during a high 
overload. This definite relationship between applied stress levels and pre- 
diction accuracy would suggest the need for a variable constraint factor which 
would be a function of the ratio for plastic zone diameter to the material 
thickness. 

Table 2 also presents crack growth predictions using no interaction 
effects, i.e., linear cumulative damage. These predictions are surprisingly 
good. This indicates that the acceleration and retardation produced by the 
load interactions nearly negate each other for most of the spectra. Only 
specimens M-88 and M-89 had a sufficiently high overload such that retarda- 
tion dominated the life. Specimens M-93 and M-94 were dominated by crack 
growth acceleration as is obvious from the very long life predicted with no 
load interaction. 
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CONCLUDING REMARKS 


This paper presented a phenomenological load-interaction model developed 
to predict crack growth propagation under any arbitrary load sequence. The 
present model is referred to as the Multi-Parameter Yield Zone (MPYZ) Load 
Interaction Model. Crack growth retardation and acceleration are accounted 
for by decreasing or increasing, respectively, the effective stress ratio 
used in a modified Forman's crack growth equation. 

The presented MPYZ load interaction model in conjunction with the 
CGR-LaRC computer routine did a good job of predicting spectrum fatigue crack 
growth in 2219-T851 aluminum. However, the predictions were not perfect as 
illustrated by M-88 and M-89 where too much crack growth retardation was pre- 
dicted. It is difficult to truly assess the prediction accurately since the 
validity of individual tests may be suspect due to a lack of replicate tests. 
Furthermore, many of the spectra were not a true test of a load interaction 
model since linear cumulative damage models worked well. 

The MPYZ model has undergone many refinements during the past 4 years 
since phenomena such as crack growth acceleration and underload effects have 
been recognized and investigated. This, indeed, increases the general appli- 
cability of this model over those models capable of accounting only for 
retardation effects. 
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TABLE 1.- MATERIAL PROPERTIES FOR 2219-T851 USED FOR PREDICTION 


CRACK GROWTH CONSTANTS 

MPYZ MODEL PARAMETERS 

C = It. 626 x 10~ 9 

A = 1.0 

n = 3.171 

B = 2.3 

K = 88 MN-m' 3/2 

Y = 0.0 

c 


AK^n = 3.3 MN/m“ 3/2 

z = 0.5 

Yield stress = 3^5 MPa 


Ultimate stress — ^tU8 MPa 





TABLE 2.- CRACK LIFE PREDICTION RESULTS 


SPECTRUM 






PREDICTIONS* 

SCALE 

SPECIMEN 






TYPE 

FACTOR 

a. 

1 

a f 

N 

TEST 

LOAD 

NO LOAD 






INTERACTION 

INTERACTION 







N /N 

PRED' TEST 

N /N 

PRED' TEST 

Air-to-Air 

0.2 

M-81 

0.16 

0.5125 

115,700 

1.08 

1.27 

Fighter 

0.3 

M-82 

0.15 

1.395 

58,585 

0.78 

0.85 

0.4 

M-83 

0.15 

0.9175 

18,612 

0.85 

0.89 

Air-to-Ground 

0.2 

M-84 

0.1575 

2.2 

268,908 

1.01 

1.09 

Fighter 

0.3 

M -85 

0.144 

1.735 

95,642 

0.78 

0.81 

0.4 

M-86 

0.1525 

1.29 

36,397 

0.71 

0.70 

Instrumentation 


M-88 

M -89 


1.805 

1.5125 

380,443 

164,738 



Navigation 

Fighter 

O C 
-p - UJ 

0.15 

0.15 

1.50 

1.40 

0.73 

0.62 

Composite 

0.2 


0.1525 


218,151 

0.99 

1.09 

Mission 

0.3 


0.15 


65,627 

0.84 

0.87 

Fighter 

0.4 


0.15 

1 

22,187 

O .90 

0.88 

Composite 

Mission 

0.2 


0.25 

0.5375 

1,354,024 

0.94 

3.00 

Transport 

0.3 


0.2575 

0.375 

279,000 

1 

0. 86 

1.8l 


Plane strain conditions assumed. 


Mean 
Std Dev 
Std Error 


0.97 

0.2k 

0.18 


1.12 

0.64 

0.38 
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cm/CYCLE 



Figure 1,- Constant amplitude crack growth rate correlation, 
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Figure 2,- Flow chart of the MPYZ load interaction model, 


27 










Figure 3,- Plastic zones at crack tip, 
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Figure 5,- Schematic of retardation in the MPYZ load 

interaction model. 




NUMBER OF TEST BLOCKS 


Figure 6,- The effect on crack growth of truncating negative 
loads from a flight-by-fl lght fighter spectrum. 




0 


Figure 7,- Schematic of underload effects in the 
MPYZ load interaction model. 
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Figure 9,- Example crack growth predictions for the fighter 

composite mission spectra, 
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